Location: Carbon Dioxide pH/Gas Transport and Exchange @ afeeca156d18 / WagnerCano SS Respiration-driven 2020-01-19.cellml

Author:
aram148 <42922407+aram148@users.noreply.github.com>
Date:
2022-07-19 16:08:46+12:00
Desc:
Added CO2 exchange and pH bg models
Permanent Source URI:
https://models.physiomeproject.org/workspace/8c0/rawfile/afeeca156d18d9435303b544ea071f7cf868595d/WagnerCano SS Respiration-driven 2020-01-19.cellml

<?xml version='1.0' encoding='UTF-8'?>
<model name="my_model" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#">
    <!-- Wagner-Cano SS respiratory model -->
    <units name="dim">
        <unit units="dimensionless"/>
    </units>
    <units name="per_mmol">
        <unit exponent="-1" prefix="milli" units="mole"/>
    </units>
    <units name="per_J">
        <unit exponent="-1" units="joule"/>
    </units>
    <units name="L_per_s">
        <unit units="litre"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="L_per_mmol">
        <unit units="litre"/>
        <unit exponent="-1" prefix="milli" units="mole"/>
    </units>
    <units name="mmol_per_s">
        <unit prefix="milli" units="mole"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="mmol_per_L">
        <unit prefix="milli" units="mole"/>
        <unit exponent="-1" units="liter"/>
    </units>
    <units name="kPa">
        <unit units="joule"/>
        <unit exponent="-1" units="liter"/>
    </units>
    <units name="mmol_per_s_per_kPa">
        <unit prefix="milli" units="mole"/>
        <unit exponent="-1" units="second"/>
        <unit exponent="-1" units="kPa"/>
    </units>
    <component name="main">
        <variable initial_value="101.33" name="P_B" units="kPa"/>
        <variable name="p_l_O2" units="kPa"/>
        <variable name="p_lc_O2" units="kPa"/>
        <variable name="p_mc_O2" units="kPa"/>
        <variable name="p_m_O2" units="kPa"/>
        <variable initial_value="0.03333" name="p_m_50" units="kPa"/>
        <variable initial_value="3.573" name="p_50_O2" units="kPa"/>
        <variable initial_value="0" name="V_In_air" units="L_per_s"/>
        <variable initial_value="1" name="V_In_air_max" units="L_per_s"/>
        <variable initial_value="2.842" name="Q_m_O2max" units="mmol_per_s"/>
        <variable name="Q_l_O2" units="mmol_per_s"/>
        <variable name="Q_m_O2" units="mmol_per_s"/>
        <variable name="C_lc_O2" units="mmol_per_L"/>
        <variable name="C_mc_O2" units="mmol_per_L"/>
        <variable initial_value="8.615" name="C_O2_max" units="mmol_per_L"/>
        <variable initial_value="0" name="t" units="dim"/>
        <variable initial_value="0.2093" name="F_In_O2" units="dim"/>
        <variable initial_value="3" name="n" units="dim"/>
        <variable initial_value="22.4e-3" name="f_O2" units="L_per_mmol"/>
        <variable initial_value="0.383" name="V_b" units="L_per_s"/>
        <variable initial_value="0.2846" name="D_l_O2" units="mmol_per_s_per_kPa"/>
        <variable initial_value="0.5692" name="D_m_O2" units="mmol_per_s_per_kPa"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>V_In_air</ci>
                </apply>
                <apply>
                    <times/>
                    <cn cellml:units="dim">1.0</cn>
                    <ci>V_In_air_max</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Q_l_O2</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>V_In_air</ci>
                        <apply>
                            <minus/>
                            <ci>F_In_O2</ci>
                            <apply>
                                <divide/>
                                <ci>p_l_O2</ci>
                                <ci>P_B</ci>
                            </apply>
                        </apply>
                    </apply>
                    <ci>f_O2</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Q_l_O2</ci>
                <apply>
                    <times/>
                    <ci>D_l_O2</ci>
                    <apply>
                        <minus/>
                        <ci>p_l_O2</ci>
                        <ci>p_lc_O2</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Q_l_O2</ci>
                <apply>
                    <times/>
                    <ci>V_b</ci>
                    <apply>
                        <minus/>
                        <ci>C_lc_O2</ci>
                        <ci>C_mc_O2</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Q_l_O2</ci>
                <ci>Q_m_O2</ci>
            </apply>
            <apply>
                <eq/>
                <ci>C_lc_O2</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>C_O2_max</ci>
                        <apply>
                            <power/>
                            <ci>p_lc_O2</ci>
                            <ci>n</ci>
                        </apply>
                    </apply>
                    <apply>
                        <plus/>
                        <apply>
                            <power/>
                            <ci>p_50_O2</ci>
                            <ci>n</ci>
                        </apply>
                        <apply>
                            <power/>
                            <ci>p_lc_O2</ci>
                            <ci>n</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>C_mc_O2</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>C_O2_max</ci>
                        <apply>
                            <power/>
                            <ci>p_mc_O2</ci>
                            <ci>n</ci>
                        </apply>
                    </apply>
                    <apply>
                        <plus/>
                        <apply>
                            <power/>
                            <ci>p_50_O2</ci>
                            <ci>n</ci>
                        </apply>
                        <apply>
                            <power/>
                            <ci>p_mc_O2</ci>
                            <ci>n</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Q_m_O2</ci>
                <apply>
                    <times/>
                    <ci>D_m_O2</ci>
                    <apply>
                        <minus/>
                        <ci>p_mc_O2</ci>
                        <ci>p_m_O2</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Q_m_O2</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>Q_m_O2max</ci>
                        <ci>p_m_O2</ci>
                    </apply>
                    <apply>
                        <plus/>
                        <ci>p_m_50</ci>
                        <ci>p_m_O2</ci>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
</model>